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A  STRONGLY  DEGENERATE  CONVECTION-DIFFUSION  PROBLEM 
MODELING  CENTRIFUGATION  OF  FLOCCULATED  SUSPENSIONS 

RAIMUND  BURGERA  AND  KENNETH  HVISTENDAHL  KARLSEN^ 


Abstract.  We  prove  existence  and  uniqueness  of  BV  entropy  solutions  of  a  strongly 
degenerate  convection-diffusion  problem  modeling  centrifugation  of  flocculated  suspen¬ 
sions.  A  modification  of  the  generalized  upwind  method  is  employed  to  solve  the  initial¬ 
boundary  value  problem  numerically,  i.e.,  to  simulate  the  centrifugation  process. 


1.  Introduction 

We  consider  the  quasilinear  strongly  degenerate  convection-diffusion  equation 

ut  +  fiu,  =  A(u)xx  +  9{u,  x),  A{u)  :=  [  a{s)ds-,  (x,  t)  G  Qt  :=  (^i,  X2)  x  (0,  T),  (1) 

Jo 

where  we  assume  that  /  >  0,  supp/(-,x)  C  [0,1],  a[u)  =  0  for  u  <  Uc  and  u  >  1,  and 
a(u)  >  0  otherwise,  i.e.,  equation  (1)  is  of  hyperbolic  type  for  u  <  Uc  and  u  >  1  and  of  para¬ 
bolic  type  for  Uc  <  u  <  1.  We  assume  that  /(-,  x)  is  continuous  and  piecewise  differentiable 

with  ||5u/(-,2;)||  <  M  and  that  A(-)  and  ^(^x)  are  Lipschitz  continuous  uniformly  in  x, 

and  that  f{u,  ■),g{u,  ■)  G  C^(xi,X2)  uniformly  in  u.  In  particular,  the  diffusion  coefficient 
a(-)  is  allowed  to  be  discontinuous.  We  consider  the  initial  and  boundary  conditions 

u(x,  0)  =  Uo(x),  0  <  Uo(x)  <  1,  X  G  [xi,X2],  (2) 

/(u,  Xb)  -  A(u(xb,  t))x  =  0,  Xb  G  (xi, X2},  t  G  (0,  T],  (3) 

and  assume  that  the  initial  function  satisfies 
uo  G  {u  G  BV{xi,X2)  :  u{x)  G  [0,1];  3Mo  >  0  :  Ve  >  0  :  IV {^xux2){^xA^ [u))  <  Mq},  (4) 
where  A^  is  defined  in  terms  of  a  standard  C°°  mollifier  with  suppcug  C  {-£,£)  via 

a^(u)  :=  ((o -f  e)  *  a;e)(u),  A^{u)  :=  (  a^{s)ds.  (5) 

Jo 

Remark  1.  If  A  e  C^,  then  it  is  sufficient  to  assume  that  TV(a:i,a;2)('^o)  <  Moreover, 
the  regularity  assumption  on  f{u,  •)  used  this  paper  can  be  weakened,  see  [14]  for  details. 
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The  Study  of  strongly  degenerate  hyperbolic-parabolic  equations  such  as  (1)  is  in  part 
motivated  by  a  recent  theory  of  sedimentation-consolidation  processes  [8].  While  the  paper 
[4]  is  concerned  with  the  case  of  settling  in  a  gravitational  field,  which  leads  to  an  equation 
similar  to  (1)  but  without  a  source  term,  we  here  focus  on  the  application  of  this  theory  to 
a  centrifugal  field  in  a  rotating  frame  of  reference  [3].  It  is  the  purpose  of  this  contribution 
to  briefly  outline  how  the  analysis  of  strongly  degenerate  convection-diffusion  problems 
presented  in  detail  in  [4,  14]  (see  also  [6])  can  be  extended  to  the  initial-boundary  value 
problem  (IB  VP)  (l)-(3),  and  to  draw  attention  to  a  new  application  of  strongly  degenerate 
equations.  In  addition,  we  utilize  an  adaptation  of  the  generalized  upwind  finite  difference 
method  presented  in  [5]  to  solve  the  IBVP  numerically,  that  is  to  simulate  the  centrifugation 
process.  For  an  overview  of  the  mathematical  and  numerical  theory  of  (strongly)  degenerate 
parabolic  equations,  we  refer  to  the  lecture  notes  by  Espedal  and  Karlsen  [11]. 

This  paper  is  organized  as  follows.  In  Section  2  we  state  and  comment  the  definition 
of  entropy  solutions  of  the  IBVP.  We  then  outline  in  Sections  3  and  4  the  existence  and 
uniqueness  proofs  of  entropy  solutions,  following  the  vanishing  viscosity  method  and  re¬ 
cent  ideas  by  Carrillo  [9]  and  Karlsen  and  Risebro  [14],  respectively.  We  come  back  to  the 
application  to  centrifugation  in  Section  5,  in  which  we  present  a  working  numerical  algo¬ 
rithm  and  numerical  solutions  for  the  IBVP. 


2.  Entropy  solutions 

Due  to  both  the  degeneracy  of  the  diffusion  coefficient  o(-)  and  the  nonlinearity  of  /(•,  x), 
solutions  of  (1)  are  discontinuous  and  have  to  be  defined  as  entropy  solutions. 

Definition  1.  A  function  u  €  L°°{Qt)  n  BV{Qt)  is  an  entropy  solution  of  the  IBVP 
(l)-(3)  if  (a)  A{u)^  G  (b)  for  almost  all  t  €  (0,T),  7xb{/(''^u)  -  ^(«)x)  =  0, 

rcb  €  {xi,X2},  (c)  lim^o  ^^(3;,  t)  =  Uo{x)  for  almost  all  x  €  (a;i,rr2),  and  (d) 


Vif  e  C^{Qt),  eR:  JJ  ^\u-k\(pt 

+  sgn(u  -  k)  [f{u,  x)  -  f{k,  x)  -  A{u)^]  -  [f^{k,  x)  -  g{x,  u)]  (p  |  dtdx  >  0.  (6) 


Here  7^,^  and  ^X2  denote  the  traces  with  respect  to  x  i  Xi  and  x  f  X2,  respectively. 
Entropy  inequalities  such  as  (6)  go  back  to  Kruzkov  [16]  and  Vol’pert  [17]  for  first  order 
equations  and  to  Vol’pert  and  Hudjaev  [18]  for  second  order  equations. 

Remark  2.  The  BV{Qt)  assumption  in  Definition  1  is  only  used  to  ensure  the  existence 
of  the  traces  jxi  o,nd  jx2  ■  Moreover,  note  that  we  can  require  that  the  boundary  conditions 
are  satisfied  in  a  pointwise  sense  almost  everywhere,  whereas  Dirichlet  boundary  conditions 
such  as  those  stated  in  [4]  (which  are  not  considered  here)  have  to  be  treated  as  entropy 
boundary  conditions  [7].  Finally,  we  point  out  that  it  is  at  present  not  known  whether  jump 
conditions  for  hyperbolic-parabolic  equations  such  as  those  by  Wu  and  Yin  [20]  (see  also 
[7]j  are  valid  here,  since  these  jump  conditions  rely  on  stronger  regularity  properties  (for 
example,  Lipschitz  continuity)  of  the  diffusion  coefficient  a(-)  than  is  stipulated  here. 
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3.  Existence  of  entropy  solutions 


Existence  of  entropy  solutions  is  shown  by  the  vanishing  viscosity  method.  To  this  end, 
consider  the  regularized  uniformly  parabolic  IB  VP,  in  which  the  functions  /,  g  and  uq 
have  been  replaced  by  smooth  approximations  that  ensure  compatibility  conditions  and 
existence  of  smooth  solutions  for  fixed  £  >  0,  and  where  is  defined  in  (5): 

+ riu,x)x  =  A^{u^)xx  +  9"iu\x),  {x,t)eQT,  (7a) 

u^{x,0)  =  ul{x),  X  e  {xuX2y,  /^(u^Xb)  -  A^{u^{xb,t))x  =  0,  x^e  {xi,X2},  t  e  (0,T]. 

(7b) 

Lemma  1.  The  following  uniform  estimates  are  valid  for  solutions  of  the  regularized 
IBVP  (7),  where  the  constants  Mi  to  M3  are  independent  of  e: 

<  -^i;  <  M2  and  ||Wi(-, t)|lL-(xi,i2)  <  M^Vt  e  (0,T].  (8) 


Sketch  of  proof.  The  application  of  the  maximum  principle  to  establish  the  first  estimate 
is  standard  and  not  repeated  here;  we  refer  to  the  proof  of  Lemma  9  in  [4]  for  details. 

Defining  sgn^(r)  :=  sgn(r)  if  |t|  >  p,  sgn^(r)  r/p  if  |r|  <  p  and  \x\r,  :=  sgn^{f)df 
for  77  >  0,  we  obtain  by  differentiating  (7a)  with  respect  to  x,  multiplying  the  result  by 
sgn^('u|),  integrating  over  0  <  Tq  <  T,  integrating  by  parts  and  using  (7b) 


/*X2  A“X2  pTo  rp 

/  \ul(^,To)\  dx  <  \{uiy{x)\  dx+  sgn^{ul{x,t))ul{x,t)  dt+  sgn(,«) 

Jxi  j  X\  JQ  d  J  Qtq 

-  dx[a^{u^))]uldtdxP  [[  sgn^{ul){gl{uyx)ul+gl{uyx))dtdx.  (9) 

dd  Qtq 


To  derive  the  second  inequality  of  Lemma  1  from  (9),  we  repeat  the  proof  of  part  (a)  of 
Lemma  11  in  [4]  to  estimate  the  first  three  integrals  of  the  right-hand  side  of  inequality  (9) 
for  77  — >■  0.  The  integrand  of  the  last  term  can  for  77  -A  0  be  rewritten  as  -|-sgn(7i^)(7^. 
Since  g^  and  g^  are  uniformly  bounded  due  to  our  assumptions  on  g,  the  desired  estimate 
on  ||w|(-,  t)||z,oo(xi,a:2)  can  be  established  by  an  application  of  Gronwall’s  lemma. 

The  same  argument  can  finally  be  employed  to  extend  the  derivation  of  the  estimate  on 
ll^t(’> ^)IU°°(xi,x2)  ici  Lemma  11  of  [4]  to  the  present  equation  with  source  term.  ■ 


From  the  estimates  established  in  Lemma  1  we  may  conclude  that  there  exists  a  sequence 
£  =  £n  iO  such  that  the  sequence  of  solutions  {u^"}  of  the  IBVP  (7)  converges  in  M{Qt) 
to  a  function  u  €  L°^[Qt)  H  BV{Qt).  We  now  have  to  show  that  u  is  actually  an  entropy 
solution  of  the  IBVP  (l)-(3).  Part  (a)  of  Definition  1  follows  from  the  following  lemma, 
whose  (short)  proof  is  a  straightforward  extension  of  that  of  Lemma  10  in  [4]: 

Lemma  2.  The  limit  function  u  of  solutions  of  (7)  satisfies  A{u)x  £  L'^[Qx). 

Finally,  repeating  the  proofs  of  Lemmas  5  and  12  of  [4],  we  can  show 
Lemma  3.  The  viscosity  limit  function  u  of  solutions  u^  of  (7)  satisfies  (6)  and  the  initial 
and  boundary  conditions  mentioned  in  Definition  1. 

Summarizing,  we  have: 
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Theorem  1.  The  IB  VP  (l)-(3)  admits  an  entropy  solution  u. 


4.  Uniqueness  of  entropy  solutions 

After  the  important  work  by  Carrillo  [9],  the  uniqueness  proof  for  entropy  solutions 
of  degenerate  parabolic  equations  has  become  very  similar  to  the  “doubling  of  variables” 
proof  introduced  by  Kruzkov  [16]  for  hyperbolic  equations.  The  key  proposition  allowing  to 
apply  Kruzkov’s  “doubling  device”  to  second  order  equations  is  the  following  version  (see 
Karlsen  and  Risebro  [14]  for  its  proof)  of  an  important  lemma  from  [9],  which  identifies  a 
certain  entropy  dissipation  term  (i.e.,  the  right-hand  side  of  (10)  below). 

Lemma  4  (Carrillo’s  lemma).  Let  u  be  an  entropy  solution  of  the  IBVP  (l)-(3).  Then, 
for  any  non-negative  if  G  C^{Qt)  and  any  k  E  (uc  1),  we  have 


jj  (^\u  -  k\ipt  4- sgn{u  -  k)  [f{u,x)  -  f{k,x)  -  A{u)^]ipx 

-  [Ixik,x)  -  g{x,u)\g)  )  dtdx  =  hm  Jj  (A(u)J^sgn(,(A(u)  -  A{k))(pdtdx. 


(10) 


Equipped  this  lemma,  one  can  prove  the  following  main  theorem: 

Theorem  2.  If  v  and  u  are  two  entropy  solutions  of  the  IBVP  (l)-(3),  then  we  have  for 
any  (p  G  C^{Qt)>  T  >  0- 

[[  (\v-  -l-sgn(u  -  u)[f{v,x)  -  f{u,x)  -  (A(r;)j,  -  A{u)y)]ipx 
J  J  Qt  (11) 

-h  sgn(r;  —  u)  \g{v,  x)  —  g{x,  u)]  dtdx  >  0. 

Sketch  of  proof.  The  argument  given  below  relies  on  Lemma  4  and  Kruzkov’s  idea  of 
doubling  the  number  of  dependent  variables  together  with  a  penalization  procedure.  We 
let  V?  €  C°°{Qt  X  Qt),  T>0,<p  =  ip{x,t,y,s),  v  -  v{x,t),  u  =  u{y,s),  and  introduce  the 
’’hyperbolic”  sets  ^  {ix,t)  E  Qt  ■  v{x,t)  <  u,,  or  v{x,t)  >1}  associated  with  v  and 
=  {(y,  s)  eQt  ■  u{y,  s)  <  Uc  or  u{y,  s)  >  l}  associated  with  u. 

From  the  entropy  inequality  for  v{x,t)  (with  k  =  u{y,t)),  the  entropy  inequality  for 
u{y,s)  (with  k  =  v{x,t)),  and  Lemma  4,  the  following  inequality  was  derived  in  [14]: 


QtxQt 


JJJJ  ~  '^1  {Pi  sgn(t;  —  m)  [/(d,  x)  —  / (u,  y)  —  (A(u)x  —  A(u)y)]  (^cpx  +  p>y) 

T 

4-  sgn(i;  -  u)  [g{v,x)  -  giujy)]^)^  dtdxdsdy  -f-  Econv 
>  lim  ////  (A(u):r  -  A(u)j/)^sgn(^(A(r;)  -  A{u))ip  dtdxdsdy  >  0,  (12) 


(Qt\£Ax(Qt\£v) 
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where  the  “error  term”  Ecnnv  takes  the  form 


Econv  =  jjjj  sgn(i;-u)(^[(/(w,y)-/(w,a:))(^]^-  [{f{v,x)-f{v,y))(p]yj  dtdxdsdy. 
QtxQt 

We  are  now  on  familiar  ground  [16]  and  introduce  in  (12)  the  test  function 

where  {5/J/i>o  is  a  standard  regularizing  sequence  in  R.  Observe  that 

Vt  +  >Pt  =  'P,  +  <P,  = 

Following  [14]  (see  also  [16]  since  f{u,  •)  is  smooth),  one  can  show  that  lim/14,0  Econv  =  0. 
Consequently,  by  sending  h  |  0  in  (12),  we  get  (11).  ■ 


Define  for  sufficiently  small  h>  0  the  functions 


Ph{x)  := 


/X 

Phix) 

•00 


I  -  ph{x  -  xi -2h),  Uh{x)  :=  ph(x  -  {x2  -  2h)),  (13) 


where  {Sh}h.>o  is  again  the  standard  regularizing  sequence  in  R.  Concerning  these  functions, 
we  have  the  following  lemma  (whose  proof  is  easy): 

Lemma  5.  Let  u  €  L^{0,T-,L°°{xi,X2)).  If  the  traces  :=  {'ju){xi,t)  and  7,^2^^  := 
{'yu){x2,t)  exist  a.e.  in  (0,T),  then  we  have  for  ip  €  C°°{Qt) 

Urn  JJ  da:(^(p{x,t){l- yh{x)-Ph{x))'ju{x,t)dtdx^  (^!p{xut)jj:,u  -  ip{x2,t)^x2U^  dt. 

We  are  now  in  a  position  to  deduce  from  (11)  the  following  uniqueness  result: 
Corollary  1.  Let  v,u  be  entropy  solutions  of  the  IBVP  (l)-(3)  with  initial  data  Vq,uo, 
respectively.  Then  for  allt  G  (0,T),  [[u(',  f)-u(-,  t)[[Li(j;j,a;2)  ^  exp(t[[5'|[Lip)  |[uo-'Uo|[Li(a;i,2;2)- 
In  particular,  the  IBVP  (l)-(3)  admits  at  most  one  entropy  solution. 

Proof  In  (11),  we  choose  (p{x,t)  =  ((1  -  ph{x)  -  Ph{x))x{t)  with  y  e  (7^(0,  T),  x  >  0, 
and  ph  and  Uh  defined  in  (13).  Note  that  ip  tends  to  xit)  as  h  i  0.  Taking  the  limit  h  |  0, 
we  obtain  from  Lemma  5  and  the  boundary  conditions  at  x  =  a,  b  (see  Definition  1): 


u  —  v\x'{t)  dtdx  > 


From  (14),  it  follows  that 


/ /  sgn(v  -  u)  [g{v,  x)  -  g{u,  a:)]  x(t)  dtdx. 
J  J  Q'p 


(14) 


u  -  v\x'{t)  dtdx  <  l|c/[|Lip 


\u  —  u]x(t)  dtdx. 


(15) 


Fixing  r  E  (0,T),  choosing  x{t)  as  ph{t)  -  ph{t  -  r)  in  (15),  subsequently  sending  h  4.  0, 
and  using  Gronwall’s  lemma,  we  get  the  L^  stability  estimate  [|u(-,r)  -  u{-,t)\\li{xi,x2)  ^ 
exp(r[[^||Lip)[|u(-,0)  -  w(-,0)[[ii(,ri,,,2)-  Since  r  E  (0,T)  was  arbitrary,  we  are  finished.  ■ 
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a)  7  =  0  b)  7  =  1 


Figure  1.  '  (a)  Rotating  tube,  (b)  rotating  axisymmetric  cylinder. 

Remark  3.  Following  [10,  13,  14],  it  is  possible  to  prove  that  the  unique  entropy  solution 
of  the  IBVP  (l)-(3)  depends  continuously  on  the  nonlinearities  in  the  problem. 

5.  Batch  centrifugation  of  flocculated  suspensions 

5.1.  Introduction.  Extending  the  arguments  developed  in  [8]  for  a  purely  gravitational 
force  to  rotating  systems,  Burger  and  Concha  [3]  show  that  the  settling  of  a  flocculated 
suspension  in  a  tube  with  constant  cross  section  (7  =  0,  Fig.  1  (a))  or  in  a  rotating 
axisymmetric  cylinder  (7  ==  1,  Fig.  1  (b))  is  governed  by  the  field  equation 

dt^P  +  r-WrifMWr^^^)  =  r-'^dr{rWrA{<j>)),  (16) 

in  which  0  is  the  sought  volumetric  solids  concentration,  r  the  radius,  the  centrifugal 
Kynch  flux  density  function,  w  the  angular  velocity  and  A(-)  is  defined  as  in  (1).  The 
function  /^.^(•)  is  a  nonnegative  Lipschitz  continuous  function  with  support  in  [0, 1],  and  the 
diffusion  coefficient  is  defined  by  a(0)  :=  (</))ag  (</>)/( Ap^),  where  Ap  >  0  is  the  solid-fluid 
mass  density  difference  and  (7g'(-)  is  the ‘derivative  of  the  effective  solid  stress  function.  We 
mention  that  equation  (1)  is  based  on  the  neglection  of  both  the  gravitational  and  Coriolis 
forces  compared  to  the  centrifugal  force  and  refer  to  [3]  for  details  on  its  derivation. 

Equation  (16)  inherits  its  degeneracy  from  the  constitutive  assumption  that  crg(-)  vanishes 
while  the  solid  floes  are  not  in  touch  with  each  other,  i.e.,  while  the  local  concentration 
(j)  does  not  exceed  a  critical  value  0c,  and  that  (7g'(0)  >  0  for  0c  <  0  <  1-  Assuming  for 
simplicity  that  supp/^^  =  (0,1),  we  see  that  a(0)  =  0  and  hence  (16)  is  of  hyperbolic 
type  for  0  <  0c  and  0  >  1  and  that  otherwise  a(0)  >  0,  i.e.,  (16)  is  of  parabolic  type. 
Consequently,  (16)  is  indeed  strongly  degenerate.  The  special  case  a  =  0  is  included  in  our 
discussion  and  corresponds  to  the  equation  studied  by  Anestis  and  Schneider  [1]. 

We  assume  that  r  varies  between  an  inner  and  outer  radii  Rq  >  0  and  R  >  Rq.  The 
solids  phase  velocity  vanishes  at  r  =  and  r  =  R,  which  implies  the  boundary  conditions 

(4(0)a;2Rb  +  5rA(0))(Rb,t)  =  O,  t>0,  R^e{R,Ro}.  (17) 

The  initial  condition  is 


0(r,  0)  =  (j)o{r),  RQ<r  <R. 


(18) 
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Differentiating  out  the  convection  and  diffusion  terms  in  equation  (16),  we  obtain  in  view 
of  the  model  assumptions  an  IBVP  of  the  type  (l)-(3).  The  existence  and  uniqueness 
analysis  therefore  states  that  the  centrifugation  model  admits  a  unique  entropy  solution  (p. 


5.2.  Numerical  algorithm.  To  solve  the  IBVP  given  by  (16),  (17)  and  (18)  numerically, 
we  employ  a  modification  of  the  generalized  upwind  finite  difference  method  presented  in 
detail  in  [5]  for  gravity  settling.  For  an  overview  of  numerical  methods  for  approximating 
entropy  solutions  of  degenerate  parabolic  equations,  we  refer  again  to  [11]. 

Let  J,  TV  G  N,  Ar  :=  {R  -  Ro)/J,  At  :=  T/N,  r,-  :=  Ro  +  jAr,  j  =  1/2, 1, 3/2, . . . ,  J  - 
1/2,  J  and  ~  (j){rj,nAt).  The  computation  starts  by  setting  :=  0o(^'j)  for  j  — 
0, . . . ,  J.  Assume  then  that  values  (/>”,  j  =  0, . . . ,  J  at  time  level  :==  nAt  are  known.  To 
compute  the  values  we  first  compute  the  extrapolated  values  :=  (p'j  —  (Ar/2)s^ 
and  :=  (p'j  +  (Ar/2)s”  for  j  =  1, . . . ,  J  -  1,  where  the  slopes  s”  can  be  calculated,  for 
example,  by  the  minmod  limiter  function  M(-,  •,  ■)  in  the  following  way: 

s”  =  )/2,  -  «ip)/Ar,  i  =  2, . . . ,  J  -  2.  (19) 

where  MM(a,  b,  c)  =  min{o,  b,  c}  if  a,b,c>  0,  MM(a,  b,  c)  =  max{a,  5,  c}  if  a,  6,  c  <  0  and 
MM(a,  b,c)  =  0  otherwise.  Moreover  we  set  Sq  =  s”  =  Sj_-^  =  s"  =  0. 

The  extrapolated  values  (p^  and  (pf  appear  as  arguments  of  the  numerical  centrifugal 
Kynch  flux  density  function  /cf°(-,-)  which,  according  to  the  Engquist-Osher  scheme,  is 
defined  by  :=  /+(m)  +  where  /+(u)  :=  4(0)  +  /“ max {4(s),  0}  ds  and 

4(i;)  ;=  min{4(-s))0}  ds.  The  interior  scheme,  which  approximates  the  field  equation 

(16)  and  from  which  the  values  (pi,.. (p'}_i  are  calculated,  can  then  be  formulated  as 


n+l  _ 

'j  — 


€ 


7-7  Ar 


(20) 


+  -  Mm  -  rUMm  -  /  =  1,  ■  ■  ■ ,  i  - 1- 


The  boundary  formulas  follow  by  considering  (20)  for  j  =  0  and  j  -  J  and  inserting  the 
discrete  versions  of  the  boundary  conditions  (17).  This  leads  to 


n  — 1 
0 


..,1+7  /EO 

RlA, 


(21) 

(22) 


To  ensure  convergence  of  the  numerical  scheme  to  the  entropy  weak  solution  of  the  IBVP, 
the  CFL  stability  condition  Rio'^  max^  14  (*?^)  I  (At/Ar)  +  2max0a((/))(At/Ar^)  <  1  must  be 
satisfied.  In  this  work,  this  condition  was  ensured  by  selecting  Ar  freely  and  determining 
At  appropriately.  The  accuracy  was  J  =  400.  For  more  details  about  the  upwind  method 
and  its  convergence  analysis,  we  refer  to  [12,  15]. 


5.3.  Numerical  example.  Sambuichi  et  al.  [19]  published  centrifugation  experiments 
with  three  different  flocculent  aqueous  suspensions,  namely  of  limestone,  yeast,  and  clay, 
using  a  cylindrical  centrifuge.  For  each  material,  the  measured  gravitational  settling  rates 
led  to  a  function  4(^)’  compression  data  determined  a  unique  effective  solid  stress 
function  o^{<p)  for  each  material.  In  this  paper,  we  choose  the  published  data  referring 
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Figure  2.  The  functions  (a)  /ck (</*))  (b)  cre{(t>)  (c)  the  resulting  function 
a(0)  used  for  simulating  centrifugation  of  a  limestone  suspension. 


to  a  limestone  suspension  (see  [3]  for  the  case  of  a  clay  suspension).  Sambuichi  et  al.  [19] 
approximated  the  measured  gravity  settling  rates  for  various  initial  concentrations  by  three 
different  connecting  straight  segments  in  a  logarithmic  plot,  which  yields  the  function 


fM)  =  < 


V 


(-47.923(/>2  +  2.5474(^)  x  10"^  m/s 
1.3580  X 

5.6319  X  m/s 

5.9735  X  10-^^-^-^^  m/s 

0 


for  0  <  (/>  <  0.035, 
for  0.035  <  (f>  <  0.08, 
for  0.08  <4><  0.119, 
for  0.119  <  0  <  <^)max 
otherwise. 


0.45, 


(23) 


It  should  be  pointed  out  that  being  cut  at  ^zi^ax,  this  function  does  not  satisfy  all  assump¬ 
tions  on  /  stated  in  Sect.  1.  However,  due  to  the  presence  of  the  diffusion  term,  solution 
values  are  bounded  away  from  ^max,  so  the  numerical  results  presented  below  would  look 
the  same  if  the  jump  of  the  function  /ck  at  ^  =  0max  had  been  smoothed  out  in  order  to 
produce  an  example  in  which  these  assumptions  are  precisely  satisfied. 

The  solid  pressure  relationship  suggested  in  [19]  can  be  converted  into  the  function 


Oe(0)  ^  0  for  (j)  <  (f)c  :=  0.28,  a,{(f))  =  0.30184(1  -  Pa  for  cj)  >  (j>c.  (24) 


The  density  difference  for  this  material  was  Ap  =  1710kg/mT  The  functions  and 
given  by  (23)  and  (24)  and  the  resulting  diffusion  coefficient  a(-)  are  plotted  in  Figure  2. 

Figure  3  shows  numerical  solutions  of  the  phenomenological  model  calculated  with  the 
functions  (23)  and  (24)  in  the  case  of  a  rotating  cylindrical  vessel  (7  =  1).  The  left  column 
of  Figure  3  shows  numerical  settling  plots,  i.e.,  diagrams  of  iso-concentration  lines  for 
selected  values  of  (f),  and  the  right  column  displays  concentration  profiles  at  selected  times. 
The  parameters  and  the  data  that  differ  in  the  three  cases  considered,  viz.  (j)o  =  0.111  and 
(jj  =  146.4  rad/s;  (?!>o  =  0.138  and  w  =  146.4  rad/s;  0o  =  0.138  and  u  =  104.9  rad/s,  were 
chosen  in  such  a  way  that  the  simulated  supernate-suspension  interfaces  could  be  compared 
with  measurements  by  Sambuichi  et  al.  [19],  which  are  shown  as  open  circles  (o).  Figure  3 
thus  illustrates  the  different  effects  of  initial  concentration  and  angular  velocity  on  the 
dynamics  of  the  centrifugation  process. 
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Figure  3.  Numerical  simulation  of  the  centrifugation  of  a  flocculated  suspension. 


While  in  the  compression  zone,  where  ^  is  valid  and  hence  (16)  is  parabolic,  the  so¬ 
lutions  are  similar  to  those  of  the  pure  gravity  case  [2,  5],  there  are  some  distinctive  features 
visible  in  the  hindered  settling  zone  <  0c)  where  (16)  is  hyperbolic,  due  to  the  rotat¬ 
ing  frame  of  reference.  Most  notably,  the  vertical  iso-concentration  lines  indicate  that  the 
concentration  of  the  bulk  suspension  is  a  (decreasing)  function  of  time,  and  the  supernate- 
suspension  interface  has  a  curved  trajectory.  These  properties  have  previously  been  found 
by  Anestis  and  Schneider  [1],  who  determined  exact  solutions  to  the  centrifugation  model 
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under  the  assumption  that  =  0,  i.e.  ^  =  0,  using  the  method  of  characteristics.  Of 
course,  in  the  centrifugal  case  (in  contrast  to  the  gravitational)  characteristics  are  not 
iso-concentration  lines,  see  [1,  3]. 
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